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In this paper the possibiUty of generating large scale curvature perturbations induced from the 
entropic perturbations during the waterfall phase transition of standard hybrid inflation model is 
studied. We show that whether or not appreciable amounts of large scale curvature perturbations 
are produced during the waterfall phase transition depend crucially on the competition between the 
classical and the quantum mechanical back-reactions to terminate inflation. If one considers only the 
classical evolution of the system we show that the highly blue-tilted entropy perturbations induce 
highly blue-tilted large scale curvature perturbations during the waterfall phase transition which 
dominate over the original adiabatic curvature perturbations. However, we show that the quantum 
back-reactions of the waterfall field inhomogeneities produced during the phase transition dominate 
f^ , completely over the classical back-reactions. The cumulative quantum back-reactions of very small 

^\) . scales tachyonic modes terminate inflation very efficiently and shut off the curvature perturbations 

evolution during the waterfall phase transition. This indicates that the standard hybrid inflation 
model is safe under large scale curvature perturbations during the waterfall phase transition. 
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I. INTRODUCTION 



Inflation proved to be very successful both theoretically [l[ and observationally [2| as a theory of early universe. The 
Oh' simplest models of inflation consist of a scalar field which is minimally coupled to gravity. A period of acceleration 
^5 ■ expansion is obtained if the potential is flat enough to allovi? for the inflaton field to slowly roll towards its minimum. 
With sufficient tunings in the parameters of the model, one can achieve 60 number of e-foldings or more to solve the 
horizon and the flatness problems of the standard cosmology. 
m ' Hybrid inflation [J, \^ is an interesting model of inflation containing two scalar fields, the inflaton fleld and the 

waterfall field. In Linde's original hybrid inflation, the energy density during inflation is dominated by the vacuum 
while the inflaton is slowly rolling. The waterfall field is very heavy compared to the Hubble expansion rate during 
inflation, H, and it quickly rolls to its instantaneous minimum. The potential has the property that once the inflaton 
^vj field reaches a critical value, (j) = (j)c, the waterfall field becomes tachyonic triggering an instability and inflation ends 

• quickly thereafter and the systems settles down into its global minimum. 
J— J . Usually it is assumed that the waterfall fleld does not play any role in curvature perturbations during inflation and 

f~^ ■ during phase transition. In this picture one basically borrows the technics and the results of single field inflationary 
^-H . models. That is, the super-horizon curvature perturbations, once they leave the Hubble radius, are frozen and remain 
ILJ ' unchanged until they re-enter the Hubble radius at a later time, such as at the time of CMB decouplings. 

Here we would like to examine this picture more closely. We would like to see if hybrid inflation is safe under large 
scale curvature perturbations during the waterfall phase transition. If one considers only the classical evolution of 
Vj ] the system, we show that during the phase transition the highly-blue tilted entropy perturbations can induce large 
blue-tilted curvature perturbations on super-horizon scales which can completely dominate over the original adiabatic 
curvature perturbations. However, we show that the quantum back-reactions of the waterfall field inhomogeneities 
produced during the phase transition become important before the classical back-reactions become relevant. We 
demonstrate that the cumulative quantum back-reactions of the short-wavelength inhomogeneities are so strong that 
they uplift the tachyonic instability of the entropy perturbations and the curvature perturbations freezes. 

Ideas similar to this line of thought, studying the amplifications of large scale curvature perturbations during 
preheating, were studied in [6|-[13 and more recently in ^20] . 

The paper is organized as follows. In section |lT] we review the basics of hybrid inflation and obtain the background 
evolutions of the inflaton and the waterfall flelds. In section IIIII the entropy perturbations and in section IIVI their 
effects on curvature perturbations are studied. In Section |V] the classical non-linear back-reactions as well as the 
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quantum mechanical back-reactions are calculated and are compared to each other. Brief conclusions and discussions 
arc followed in section IVTl 

While our work was finished the work by Lyth [2l[ appeared which has overlaps with our results. See also |22l | 
which appeared shortly after our work. 

II. HYBRID INFLATION 

Here we study the basics of hybrid inflation 0, |j] and the background field dynamics. 

A. The Potential 

The potential in standard hybrid inflation has the form 

F(</., V) = ^ (^' - ^) ' + \m^^ + \9'<t^'^' , (1) 

where (f) is the inflaton field, V is the waterfall field and A and g are dimensionless couplings. The system has a 
global minimum given by </) = and -0 — Af/vA. Inflation takes place for (j>c < (t> < (t>i where (j)i is the initial 
value of the inflaton field and 0c = ^^ / 9 is the critical value of where the waterfall field becomes instantaneously 
massless. During inflation ip is very heavy and is stuck to its instantaneous minimum 0; ~ 0. For (j) < (pc the waterfall 
becomes tachyonic triggering an instability in the system which ends inflation abruptly. Soon after phase transition, 
the systems settles down to its global minimum and inflation is followed by the (p)reheating phase. 

As in Linde's realization of hybrid inflation [3|, we consider the limit where the inflation is dominated by the 
vacuum. For this condition to hold one requires that 



A 
? 



M^ > — m^ . (2) 



To solve the flatness and the horizon problem, we assume that inflation proceeds at least for about 60 number of 
e- foldings. In the vacuum dominated limit the number of e- foldings is given by 

2^M4 /0A 

^^^A^^^^n^j' ^'^ 

where m?p = 1/G with G being the Newton's constant. We assume that 0i is few times 0c so one can basically neglect 
the logarithmic contribution above. 

To get the correct amplitude of density perturbations, one has to satisfy the COBE normalization for the curvature 
perturbations Vu ~ 2 x 10^^. The power spectrum of curvature perturbations is 

_ I287r y3 g2 j^.jw 
^~ 3^"V^^ A3m6 m4' ^' 

where the relevant quantities are calculated at the time of Hubble radius crossing (fc — aH) at 60 e-folds before the 
end of inflation. In this picture, it is assumed that the curvature perturbations are frozen once the modes of interest 
leave the Hubble radius, as have been treated in conventional analysis of hybrid inflation so far. Our main goal in 
this paper is to examine the validity of this assumption more closely. 

We are interested in the limit where the waterfall field rolls rapidly to its global minimum once the instability is 
triggered. For this to happen, the absolute value of the ip mass should be much bigger than the Hubble expansion 
rate during phase transition so 

M^ < AmTOp . (5) 

B. The Background Fields Dynamics 

Here we study the classical evolutions of background fields and tp during infiation and phase transition, see 
also [a, I23I |24| where somewhat similar analysis were carried out too. In subsection IV Bl we study the quantum 
back-reactions to the the system in the Hartree approximation. 



The equations of motion for (j) and ip are 

4> + iH^ + (m^ + g^^?)<j) = 



(6) 

(7) 



With the ass umptio n of vacuum dominated potential, the Hubble expansion rate is nearly constant during inflation, 
H = Ha = ■\/27r/3A M"^ /nip. It is more convenient to use the number of e-foldings as the clock, dN = Hq dt and the 
background fields equations are now written as 

/ „;.2 \ 

(8) 

(9) 



(10) 



and the prime denotes the differentiation with respect to the number of e-foldings. 

Using the rapid waterfall condition Eq. ([5]) combined with the expression for the total number of e-foldings Eq. ([3]) 
one obtains 
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(11) 



Similarly, one can check that a/? 3> 1. 

To simplify the notation, we take the critical point as the reference point and define n = N — Nc- We use the 
convention that at the start of inflation for (/) = c/ii, A^ = 0, at the time of phase transition N = Nc and at the end of 
inflation N = Ne- With this convention n < before phase transition whereas n > afterwards. As assumed, the ip 
field is much heavier than Hq during inflation so it rapidly rolls down to ?/> ~ and one can simply solve Eq. [8] 



(n) ~ (/)cexp(-r n) 



with 



Equivalently, one also has 




Nc 



-ln(^ 
r 



(12) 



(13) 



(14) 



Let us now turn to the dynamics of Tp field. As it can be confirmed from our full numerical results we are in the 
limit where i^'^/Hq <^ /3/A so the equation for ip simplifies to 



V-" + 3^^' + /3 (e- 



-2r n 



i)V' = o. 



(15) 



We divide the solution into two regions. First, we solve this equation for the period before phase transition, (p > 
The solution of this equation during this period is 
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(16) 



where J^ and Yy are the Bessel functions, c\ and C2 are constants of integrations and v = ^J pi + 9/A/r. As explained 
previously, we are in the vacuum dominated limit so /? 3> 1 and u ~ ^f^/r — \f^Nf, ^ 1. Using the approximations 
for the Bessel functions with large arguments one finds 
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FIG. 1: The background classical dynamics for hybrid inflation potential Eq. ((!]) with M — .67 x 10"'', m = 2.5 x 10~^° and 
X — g — 2.5 X 10~ ^ in the units where rrip = 1. The upper solid blue curve shows the amplitude of curvature perturbations, the 
middle solid green curve shows the background (f) field evolution and the lower solid red curve shows the behavior of background 
^ field, all three curves are obtained from full numerical analysis. The dashed dark blue curve shows our analytical solution 
for ip fleld. The vertical axis is logarithmic and the horizontal axis is the number of e-foldings. 



This solution means that ^ has oscillatory behavior with exponentially decaying frequency and exponentially decaying 
amplitude. This behavior can be seen in Fig. [1] 

Now, we solve the ip evolution for the period after phase transition till end of inflation, (j) < 0c- Since inflation ends 
in few e- folds after reaches the critical point, it is more appropriate to use the small rn approximation in Eq. (jlSp 
and 



which has the following solution 



V. 



-3n/2 



V'" + 3V''- {2l3r njiP = , 



CAiiefn 



4e 



4/3 



C2Bi ef n 



4e 
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(18) 



(19) 



in which 



e^ 



/2r/3 



:«/3, 



(20) 



and Ai(x) and Bi(x) are Airy functions of first and second kind respectively. For e^ ^ 1, which is the case in our 
analysis to satisfy the water-fall condition, the second term in the argument of Airy function can be ignored. Since 
Ai(x) is a damping function for a; > we just keep the term containing Bi(x). In a good approximation the solution 
can be read as 



^(n)c^^(7V,)Bi(6^/^n) e-^"/^. 
where ip{Nc) = tpi"^ = 0). By using the asymptotic behavior of Bi(z) for large z, z — > cx), 

^2/3z^/^ 

Bi{z) ex 



and keeping just the exponential dependence one finds that for (jj < (jjc 

1 / 3 



3 — r 

ij} ~ipi exp ( — Nc 
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where to get the final result Eg. p^ have been used. 
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FIG. 2: Here the contents of Fig. [l] are shown in the last few e-foldings after the phase transition at Nc — 56. Here we see 
that, considering only the classical fields dynamics, TZ grows rapidly at the end of inflation until the condition (I56|l is met when 
7^' ~ at A'^ ~ 63.5. This is followed by a short period of inflation when condition (|57p is met and the non-linear corrections 
to '4' dynamics settles it down to the global minimum ending inflation. 



III. THE ENTROPY PERTURBATIONS 

Here we look into entropy perturbation following f25l] closely. In the field space of (0, iji) one can perform local 
fields rotations such that 



5a — cos 9 5 (j) + sin 9 dtp , Ss = — sin 6 6(p + cos 9 Sip , 



(24) 



such that cos 6' = (j) / \J (jP' + tp"^ and sin0 — ip / d (jy^ + tp"^ . In this picture, 5a and 5s represent, respectively, the 
adiabatic and the entropic perturbations. One can check that the evolution of the curvature perturbation for a mode 
with momentum k is 



where 



n = -. — ^* + 95s. 

H a^ a 



9 = -Vs/a 



(25) 



(26) 



with Vs — cos 9V^ — sin 9 V^. Eq. (j25p indicates that in the presence of large entropy perturbations or sharp turns in 
field space, the curvature perturbations can change on super-horizon scales. 
The equation of entropy perturbations is [25| 



5s + 3H5s +[ ^ + Vss+39'^\ 5s 



a 2TrGa^ 



r*. 



(27) 



in which Vss = (sin^ 9)'^V^^ - (sin26')y^^ + (cos^ 9)'^V^,^. 

We are interested in solving this equation for 9',9 -^ 1. As we shall see later, these conditions are true until the 
end of infiation which can also be seen from Fig. [3] For the later references, from Eq. (fT2)) and Eq. ((23)) . one can 
easily find that 
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/i\ exp(^n3/2-l^^ 



which for n > Q ov (j) < (pc (after phase transition) results in: 
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FIG. 3: Dynamics of angle in the {(j), ip) field space with the same parameters as in Fig. [T] : The upper solid black curve and 
the lower solid blue curve, respectively, show the full numerical solutions of |lnSj and |lnS'|. The upper dashed green curve 
and the lower dashed red curve, respectively, show the evolution of the corresponding quantities obtained from our analytical 
solutions, Eqs. (|28|l and (I29|l . Inflation ends when ~ 1 indicated in the graph by the vertical line. 



Neglecting in Eq. (P7)) . the equation governing the evolution of entropy fluctuations is the same as the equation of 
background ijj field, Eq. ([TC|) . This indicates that, although the entropy perturbations are suppressed before the phase 
transition, but they becomes highly tachyonic during the phase transition. The subsequent tachyonic enhancement of 
the entropy perturbations overcome their initial suppression prior to phase transition. This is the key effect to obtain 
large curvature perturbations induced from the entropy perturbations at the classical level. However, as we shall see 
in lVBl one should also take into account the quantum back- reactions which can uplift the tachyonic instability of the 
background "0 fields and modify the evolution significantly. 

The differences between 5s and ^ evolutions therefore are only in the initial conditions. For a given mode, we 
consider the evolution of the entropy fluctuation when it leaves the Hubble radius at the moment of Hubble radius 
crossing A^,, so for N > N^., similar to Eq. (|23l) . one obtains 



6s{n) ~ 5s^ — ^ 



1/2 



exp 



3-r 



-iV, 



exp(--7V,)exp ( - Y + ^"^^ 



(30) 



We note that the contribution exp (2e^n'^"/3) represents the tachyonic enhancement of the entropy perturbations 
during the waterfall phase transition whereas the factors exp {—3{Nc + n)/2) indicates the suppression of the entropy 
perturbations from the time of horizon crossing till the onset of phase transition. Both of these competitive behaviors 
can be seen in Fig. [T] 

To find the final amplitude of entropy perturbations, we need to find their amplitude at the time of Hubble radius 
crossing, Ss,.. To do this we note that the entropy perturbations, which are basically the waterfall field perturbations, 
are very heavy during inflation. As explained previously, this causes their suppression before the phase transition 
which should be taken into account. During inflation <C 1 so (5s = Sip. Rewriting the mode equation in conformal 
time 77 = —1/aH and defining v = a{r])ds the equation of entropy perturbation is 
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(31) 



in which 



w|(0) 



/3 K7-I 



(32) 



Since we are interested in momenta which exit the horizon during first few e-folds of inflation the above ratio is nearly 
constant and 
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(33) 



By these considerations and noting that /3 ^ 1 one has 



5V 
9772 



+ 



^2 



Wfc = 0. 



(34) 



As the frequency of this equation changes adiabaticaUy, one can use the WKB approximation with the Bunch-Davis 
vacuum for initial times or r/ — > —00 and obtains 
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Vk 



(/4fc2 + 4/3/r/2 
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(35) 



Using this solution the amplitude of the entropy perturbation at the time of Hubble radius crossing, k\r]\ = 1, is 
obtained to be 



H 1 



/2F/31/4 



(36) 



The extra factor f3 ^'^ represents the suppression of the entropy perturbations before the phase transition. 

Assuming that the scale factor at the start of inflation, N,, e-folds before the end of inflation, is unity, a{Ne) — 1, 

one also has N^ — ln(a,) == In ( — I . The final amplitude of entropy perturbations at the end of inflation therefore 



is 
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(37) 



where Uf = Ne — Nc denotes the number of e-foldings from the start of phase transition till end of inflation. For a 
vacuum dominated potential with a quick phase transition, we have n/ ^ 1. 
The final amplitude of the entropy perturbation 5 = -S-Js is 



Vs = As[ — 



(3-r) 



with 



A, 






-3iVc 



exp I — 3n/ + 



4e, 



V- 3/2 



(38) 



(39) 



From Eq. pSI ) one observes that the entropy perturbations are highly blue-tilted with n^ ~ 4 — r ~ 4. In next section 
we show that at the classical level these highly blue-tilted entropy perturbation induce large blue-tilted spectrum on 
super-horizon curvature perturbations. 

IV. POWER SPECTRUM OF CURVATURE PERTURBATIONS 

We now have all the materials to calculate the final power spectrum of curvature perturbations. For this purpose 
we need to know the amplitude of adiabatic curvature perturbations at horizon crossing as the initial conditions and 
integrate the evolutions of curvature perturbation from the time of horizon crossing till end of infiation. The final 
amplitude of curvature perturbation, therefore, is 



7^/ - 7^o 



n'dn. 



(40) 



where TZq represents the adiabatic curvature perturbations in the absence of entropy perturbations. 

Starting with Eq. (P5|). the evolution of curvature perturbations for the super-horizon modes, induced by the 
entropy perturbations, can be written as 



a' 



(41) 
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As one can see from above equation both 0' and 5s can source the curvature perturbations. We also note that 9' 
represents the acceleration of the field "0, specially during the phase transition. As can be seen from the full numerical 
analysis, our classical background is such that during inflation and phase transition, 9,9' <C 1. However, shortly 
after phase transition ip rises quickly from its value during infiation ?/; ~ to its final value at the global minimum 
tp = M/vA when inflation ends. The rapid rise of ijj and ij]' cause inflation to end when 9, 9' become large which can 
also be seen in bf Fig. [3] Equivalently, this can be interpreted as when the classical as well as quantum back-reactions 
from g'^(j)'^ip'^ and Xip^ interactions induce large masses for and i/j such that they roll rapidly to the global minimum, 
ending inflation. Therefore, in the analysis below we work in the limit where 0, 0' <C 1 and consider the end of inflation 
when 9 = 9f :^1. 

To calculate the evolution of curvature perturbation from Eq. (|^T|) we need to estimate 9' and 5s. The derivative 
of 9 in field space is 

Since r ~ 1/A"e ^ 1, the first term is much larger than the second term by a factor of e^/r and 9' ~ V"/"^'- To 
calculate 5s we observe that its equation has the same form as the background "0 equation and therefore 

5s{N) = ris^^iN) , (43) 

in which 

3-r 

5s{N^) _ (5s, f k\ 2 ^^s^ 



The term e"^*^ represents the phase difference between the oscillations of ip and the entropy perturbations. This phase 
difference vanishes after time averaging when we find the final curvature power spectrum. 

Combining the above expressions for 9' and 5s, the final curvature perturbation is integrated to 

7^(iVe) = 7^o - 2 p ^^"V dn . (45) 

As the function Bi grows more rapidly than the linear exponential, from Eq. [21] one obtains 

V^"(n) ~ i>iN,) e-3/2" Bi" (e^/^ n) . (46) 

Using Eq. (jAlOp and noting that cj)' ~ — r^c cxp(— r n) from Eq. (fT2|) . the final curvature perturbations is obtained 
to be 



7^ = 7^o - 2i;{N,f ^s^^ / " e-^^-s-)" Bi^ (e^/^^) dn . (47) 
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'r^(f>l JO 
The amphtude of quantum fluctuations of adiabatic perturbations at the moment of horizon crossing is 

TJ 



/2P 



(48) 



where Q^ represents the Sasaki-Mukhanov variables for the adiabatic perturbations [2^ . Using the form of flg^ given 
in Eq. pij) . the curvature perturbation calculated from Eq. (HTl) is 

( 3-3r \ 

1-C(^-|) ' I(n/)j , (49) 

where 7?.o is the initial value of adiabatic curvature perturbations at the moment of horizon crossing 

(50) 
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and 

r/31/4 4 

Furthermore, the integral 1{nf) has the followmg form 

I{nf) = n v! e-(3-2r)n' 3^2 |^^2^/3 ^,^ ^^ ^52) 

Eq. (j49|) has some interesting features. If the second term in the big bracket in Eq. (|49| is larger than unity, then the 
induced curvature perturbations from the entropy perturbations dominate over the adiabatic curvature perturbations. 
Furthermore, the dominant momentum dependence in the big bracket comes from the initial amplitudes of entropy 
perturbations and the integral T-inf) is nearly constant for all momenta. This implies that the induced curvature 
perturbations from the entropy perturbations are highly blue-tilted. 

Using the integral approximation given by Eq. ( IA9|) one can calculate the integral approximately and 

Z(n,)^^e-(3--)"/Bi^(efn). (53) 

Plugging this into Eq. (|49p. and using the asymptotic behavior of Airy function of second kind, the curvature 
perturbations at the end of inflation is calculated to be 



TKjif) = 7^o 




(-3+3-)^= exp (4/3e^nf 2 - (3 - 2r)n^ 



(54) 



To get the final curvature power spectrum we need to know rt/, the time of end of inflation. To determine this we 
proceed as follows. The exponential growth of the ■0 modes (background as well as quantum fluctuations) can violate 
our background solutions for (^ and '0- The exponential growth of the background V' field have two important effects. 
First, through the interaction term (7^0^-0^/2, it can induce large mass for field which speeds up its rolling toward the 
global minimum and violate the slow-roll conditions. The second important effect is that the back-reaction from A'0^/4 
term induces large positive mass for f/; which uplifts the tachyonic mass of "0 field leading to the deceleration of this 
field. These two effects jointly terminate both inflation and the growth of the super-horizon curvature perturbations. 
To see the latter effect we start from Eqs. (|42)) and (|4T|) where 



^' - ^ ■ (55) 

We observe that the deceleration of the field as well as the fast-rolling of the </> field jointly cause the termination 
of the super-horizon curvature perturbations as explained above. 

V. CLASSICAL VS. QUANTUM MECHANICAL BACK-REACTIONS 

In this section we study the back-reactions of classical (zero momentum) mode of field as well as the quantum 
back-reactions of "0 inhomogeneities produced during phase transition on the dynamics of the system. We examine 
which of the above two mechanisms dominate sooner to terminate inflation. 

We will demonstrate that the quantum mechanical back-reactions dominate completely over the classical back- 
reactions and hence the end-point of inflation is determined by quantum back-reactions. However, as we shall see 
later, the variance of the quantum fluctuations (on all scales) after phase transition has the same time dependence 
(n-dependence) as the classical trajectory and the difference is just in a proportionality factor. In order to demonstrate 
why the classical back-reactions are secondary we start the analysis with the classical back-reactions. 

A. Classical Back-reactions 

Due to smallness of background -0 field before phase transition our analysis in IIIBI concentrated only on the linear 
level. Now we add the back-reactions of the non-linear terms g'^(\P"\\P' 12 and A'0^/4 on the dynamics of the system at 
the classical level. 
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The correction from the interaction term g^cjy^ijp' /2 becomes important in </> evolution, Eq.®, when 



9 



HI 



(56) 



On the other hand, the non-hnear self-interaction term in background tp evolution, Eq. 0, becomes important when 



/3r~ A- 



7^2 



(57) 



However, Xa/g^jir ~ A//?^^ ^ 1 as a consequence of the vacuum domination condition Ea. dTTT) . This indicates 
that the back-reaction of the ip field on the inflaton dynamics becomes important sooner before its self-interaction 
corrections affect its own dynamics. Putting it another way, during hybrid inflation e <C ?7 where e and 77 are the 
conventional slow-roll parameters. Due to back-reactions of "0 on inflaton mass the condition \rj\ ~ 1 is met sooner 
before e find the chance to become order of unity due to its initial smallness during inflation. We have numerically 
checked that when the condition Eq. ((56| is satisfied then i?' ~ 0. This followed by a very short period of inflation 
when condition Eq. (j57l) is met and inflation ends. 

Combining Ea. ([55| and Eg. ip^ one finds the end-point of inflation to be 



exp 
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,3""^^ ""V ~^^ \PJ V^». 
Using this in Eq. (I54p . the amplitude of curvature perturbation at the end of inflation is 
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(59) 



Consequently, the final power spectrum of curvature perturbations, VuS^O^ ~ k') = ^p-(7?,(k)7?.(k')), is calculate to 
be 



'Pn = 'Ptzv. 



(3-3r) 



4"' 



27^2^2^2^1/2 




(60) 



where Vtzq = k^TZ^/2'r:'^ is the power spectrum of the adiabatic curvature perturbations. As explained previously, 

during inflation tp is very heavy so ip rolls to its instantaneous minimum very quickly and ipi ~ 0. Below we 

demonstrate that (pi/ipi 3> 1 such that the curvature perturbations induced from the entropy perturbations dominate 

completely over the adiabatic curvature perturbations. To see this note that for the ip field to be nearly zero such 

that its fluctuations do not contribute to the curvature perturbation around N^ e-foldings before the end of infiation, 

7^' 
that is -r^ ^ 1 at the start of infiation, one requires that 
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Using Eq. p3)) and Eq.® one can easily find that 



This in turn results in 
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Now compare the second term in the big bracket in Eq. (j60p to the initial curvature perturbations 



Vno 



h)^ 1 27r2r2/32^i/2 



4"' 




(61) 
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FIG. 4: Here the final power spectrum of curvature perturbations, at the classical level, as a fiinction of comoving momentum 
for the parameters used in Fig. [T] are shown. The momenta exited approximately in first 6 e-folds. Blue squares shows the 
full numerical results whereas the solid red curve shows the best fit. The best power law fit is acquired for V-r. = Ak"^"^^ with 
n, ~ 4.01(3.98 - 4.04) at 95% CL. 

Noting that (k/H)^ ~ e^' with N^, > 1 for super-horizon modes, and using the inequality Eq. (j63p and r ~ a/S one 
obtains 



AP„ ,^ (,jf^ir^ „_ 



T-'na 



71-2 p 



(65) 



Noting that A^* > 1 and (3 ^ (5, one concludes that ^'^ ^ N^. This indicates that, in the limit where only the 
classical back-reactions are considered to determine the endpoint of inflation, the curvature perturbations induced 
from the entropy perturbations dominate completely over the original adiabatic curvature perturbations. We shall 
see in next subsection that this conclusion is not stable against quantum back-reactions. 
With -^^ '3> 1, the spectral index of curvature perturbations, n-jz — 1 = dlnV-jz/dhik, is 



n-ji — 1 = n% — 1 + 3 — 3r , 



(66) 



in which n^ is the spectral index calculated from the adiabatic curvature perturbations using Viz , n^ — 1 = 2ri—6eci 
2ri ~ 2a/3, where e and rj are the standard slow-roll parameters. Using the approximation r ~ a/3 ~ 1/Ne one has 



"^-^^^-iv: 



(67) 



This indicates that the curvature perturbations receive a large blue-tilted spectrum from the entropy perturbations 
and n-n — i- 1/Ne. 

Finally, to calculate the time of end of inflation, Uf, from Eq. ((58)) we obtain 



Uf 



9 

4e,; 



-Nr, 



^ \n el/'C" 



4:6,1 



PJ VV^: 



-1 2/3 



(68) 



For {(l)i/'ipi) ~ lO'^/3/r, the second term in the big bracket above can be ignored and one approximately has 



"/ 



4e, 



-N, 



^ 



2/3 



(69) 



It is also interesting to calculate the angle in phase space at the end of inflation. Of. Using Eqs. 
obtain 



and 



9f ~ 3e"i''^= 



1. 



(70) 



This verify our previous claim that inflation ends classically when 9f '^ 1. This can also be seen from Fig [31 
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As an example consider the vacuum dominated hybrid inflation with parameters M = 0.67 x 10~'^TOp, m = 
2.5 X lO^^^JTip and X = g^ = 2.5 x 10^^^. From Eq. ([ST]) the spectral index is n-]^ ~ 3.99 which is in good agreement 
with n-jz — 4.01 obtained from the full numerical analysis. Also from Fig. |4]one can see that there are good agreements 
between our analytical results and the results obtained from the full numerical analysis at the classical level. Finally, 
n/ ~ 7 which is in agreement with the numerical results (see Fig. ^ Nc — bQ and A^e — 63.5 so n/ ~ 7). 

B. Quantum Fluctuations Back-reactions 

Our analysis so far concentrated on the classical evolution of the system. Due to tachyonic instability during 
the phase transition quanta of 4' particles inhomogeneities, 5V'(x, t), are produce [2y, [23] which can back-react on 
the classical backgrounds as in preheating models |28l430ll. If the mechanism of particle creation due to tachyonic 
instability is very efficient, the back- reaction of the produced particles can induce large mass for inflaton field in the 
form of g'^{Sip'^) where {Sip'^) is the expectation value of S^p in the Hartree approximation. This violates the slow-roll 
conditions and (j) rapidly rolls towards the global minimum. Furthermore, the large back-reactions of Sip also induce 
a large mass for the background ip which can uplift its tachyonic mass. Therefore, one has to take into account the 
quantum back-reactions and see which of the classical or quantum back-reactions dominate first to terminate inflation. 

To handle the quantum back-reactions, we use the Hartree approximation and calculate the effects of all modes 
which become tachyonic during the phase transition. Let us compute which modes become tachyonic after the phase 
transition. The equation of the Sipk fluctuations at the linear order is 



S^'^ + 35Vfe + ( ^j^ - 2l3rn ] Si^k ^ (71) 



Therefore modes which satisfy the inequality 



- < ^/Wr = e^ (72) 



become tachyonic soon after the phase transition. Here we defined k^ = a{Nc)Ho = e^'=Ho as the critical mode which 
leaves the horizon at the time of waterfall phase transition. 

We divide the tachyonic modes into two categories: first, large modes ipj^, corresponding to k < kc, which exit the 
Hubble radius sometime before the phase transition and second, the small modes ip^ , corresponding to fcc < A: < e^kc, 
which do not exit the Hubble radius till time of phase transition. Using Eq. (|35p , the amplitude of both modes at 
the time of critical point (j) = (pc ^s 



and therefore for large modes one has 



e-3Afc/2 



\S^Pj;{n = 0)\c^—=^. (74) 



This is in agreement with our previous result, Ea.p7p. for n = 0. On the other hand, for the small modes which 
remain sub-horizon till the time of phase transition one has 

|5V'f(n = 0)|^^=. (75) 

With this division of the modes, and using Eq. (P7)) . {SiJ;'^) is calculated to be 



(5V' 



Oft i _3^^ ; U 't ^ ^-2Na 



(27r)3 2Ho "" ' Jk^ (2^)3 2fc 



exp (-3n + ^n^/A 



= {Si;')L + {Sr)s, (76) 

in which we have ignored factors of order unity. The first integral, representing {Sip'^)^, comes from the large modes 
which leave the horizon sometimes prior to phase transition whereas the second integral, representing {Sip'^)s, indicates 
the contributions from the very small scales modes which are sub-horizon till phase transition. 
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We observe that both {Sip'^)s and {SiP'^)l have similar n-dependence which becomes important in our discussion 
below when we compare the classical and quantum back-reactions. However, 

(<5V/)s^4W)l- (77) 

As e5, ^ 1 the second integral in Eq. ([7^ is much bigger than the first one. This means that the cumulative 
contributions of the modes which becomes tachyonic but remained sub-horizon till the time of phase transition are 
more important in quantum back-reactions and 

i^^')..i^i^')s_ 4 exprifin3/^-3.V (78) 



ij2 i/2 4^2 -"^\^ 3 

As explained before, there are two competitive corrections which terminate inflation and the evolution of large scale 
curvature perturbations during the phase transition. The First one is the classical corrections from the (7^0^^^/2 
interactions which induce large mass for the inflaton, violating the slow-roll condition during the phase transition. 
This is the effect which was studied in subsection (jV Al) . specifically Eq. ([55)) . The second contribution is the quantum 
back-reaction corrections into inflaton mass via g'^(j>'^ (Sip^)- We need to see which of the above two corrections dominate 
first. Comparing Eq. ()58|) and Eq. ()78|) . one finds that these two terms have similar time-dependence (n-dependence) 
and 



(S^^n)) {S^'{n))s 2 9 



Using Eq.([TT]) and Eq.®, one finds that 



4-F - e^"=- (79) 



■02 (n) 



(g2/33/2) e^^^^^ . (80) 



As exp(3Afc) ^ 10^^, one concludes that for any reasonable value of coupling g the quantum back- reactions dominate 
completely over the classical back-reactions to terminate inflation. Therefore, the end-point of inflation is determined 
by the quantum back- reaction effects from {Sip'^{n))s and 

4 3/2 ^ N STT^a 



exp(-e^n/ -Sufj--:^. (81) 

It's now instructive to calculate Uf, the time when the quantum back-reactions terminate inflation. If we follow 
the estimation of [2^, [23 performed in flat backgrounds, we obtain 

1 , Sn'^a , , 

Hf In —^ . 82 

For the parameters of our numerical investigations such as in Fig [1] one has Uf ~ 3. However, taking into account 
the background cosmological expansion, from our Eq. (|8T|) , one flnds 



which for the parameters of our numerical studies this gives Uf ~ 1.7. This clearly indicates that quantum back- 
reactions dominate sooner than the classical back-reaction which happens at nj ~ 7. 

Another important point which should be considered is that the quantum fluctuations, ((5'0^), can change the 
effective classical trajectory. As we demonstrated in Eq. ([79]) the expectation value of the quantum fluctuations 
dominates completely over the "zero momentum" mode and they will induce new effective trajectory in the field space. 
The fact that (S-ip^) dominates over the background classical field contributions in waterfall dynamics also indicates 
that using the ip "zero momentum" as the classical trajectory is not reliable. Quantum modes with tachyonic mass 
become highly occupied soon after the phase transition and in some senses they become classical. This suggests that 
we may take ^JJEilP) as the effective classical field trajectory. Below we justify this proposal so one can introduce an 
effective classical trajectory defined by ^ \f(M^ ■ 
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To justify our proposal, we start with the conventional method for the evolution of super-horizon curvature pertur- 
bations [43II. In the study of curvature perturbations in hybrid inflation this method was pioneered in [2l| and [ij] 
(see also [42|)- In this method, the change in the comoving curvature perturbation on super-horizon scales is given by 

7^ek-i/%^, (84) 

where 5pc is the pressure perturbations on comoving slices. In our formalism we relate 5pc to the entropy perturbations 

5s or Sip. In the analysis here it is assumed that the ip field is frozen at the background so '0 = classically. Our goal 

is to demonstrate that on super-horizon scales Eq. (|84p reduces to our starting equation for the evolution of curvature 

perturbation, Eq. (j25p . obtained in [25] for the two field inflationary system, with the appropriate definition of 9 and 

replacing if) — >• y/{6i{j'^). 

One can simply check that p + p = d^, where taking into account the quantum effects, one also has 42] a^ = 
• 2 
0^ -|- ((S?/; ). However, as in our previous analysis, the velocity along the infiation trajectory in the field space can be 

well approximated by &^ — (fP' till end of infiation. Therefore 

Tte.^^^. (85) 

a a 

As the potential is vacuum dominated we can also neglect the gravitational back-reactions on Sip. Neglecting the self 
interaction term, and noting that the comoving slice coincides with the (50 = surface [42] , the contribution of ^- field 
in the pressure becomes 

Jp, = ijV^%i(Af2-gV)<5V''- (86) 

Using Eq. [3D] for the time evolution of ds mode (which is the same as 6i> mode) one finds that the contribution of 
the kinetic term above is approximately equal to the potential term and 

5p,^c^{M^-g^cj,^){5i;^)^. (87) 

Below we would like to find an expression for ((5?/^^). to calculate 5pc\s.- One can show that [42] 

Most of the contributions to the momentum p integration above comes from the small scales modes with the amplitude 

5Vk(n) = 5Mn = 0)f{n) , (89) 

in which f{n) is defined via 

f{n) ^ exp ten3/2 - Z^ , (90) 

and (5V'k("- = 0) = e^^' j\r2k as given by Eq. ([75]) . Performing the integral and setting the UV cutoff p = e^fcc yields 

= {27TfS^k + k')-^-^ J dpp\6^Pl\ fin) 

= {27rfS^{k + k')^-^e^Hof\n). (91) 

On the other hand, for the large scale quantum fiuctuations we have 

iSi^tSi^t') = {2nf6'{k + k')— — /(n)2 , (92) 

zhq 

where Eq. (|74|) have been use for their amplitudes at the time of waterfall. Furthermore, as we showed in Eq. [78| 

W') = §4/N' • (93) 
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Combining Eqs. ^, dMI and ^ we obtain 

{{S^^)^ {S^^)^,) = -{5i;'){Si;tS^t') ■ (94) 

This equation suggests that, as long as we are interested in two-point functions at the linear perturbation theory, one 
can make the following identification 




ViWJHt ■ (95) 

Having obtained this identification for ((5-0^), we plug it into Jpck expression in Eq. ([57)) and Eq. (|55|) to obtain 



T^ck = ^\ 7^ ■■ — '5'/'k • (96) 



Now we can cast Eq. ([96]) in the form of our starting formula Eq. p5|) for TZ as prescribed in[25|. Noting that 
(M^ — g'^<t>^)\/ {S^l)"^) = — y^(f/' — > \/W4^) and (5sk = Sipy^ on super-horizon scales, Eq. (|96l) can be written as 

nu = '^e,sSskJ^ (97) 

where similar to Eq. ([26)) 

fl V 



(98) 



Therefore we have justified our proposal in using y^ {6ijj'^) as the effective trajectory when the background ip field is 
zero. Compared to our starting formula Eq. ([^5t . there is an extra factor yjxjlt^, in Eq. ([W)) which originated from 
the UV cutoff imposed in Eq. ([91]) . However, this extra factor is not significant and it does not affect our analysis 



below. For instance, in our numerical example we have ^^1/26^ ~ 1/4. 

The quantum back-reactions therefore have two crucial effects. First they determine the end of inflation given by Eq. 
([79)) and Eq. ([8T)) . Second, as demonstrated above, they provide the effective classical trajectory via i/; -^ \/((5'0^). 



To see how the replacement V' — > yWiP^ changes our results in previous sections based on classical analysis we note 
that the only place in which we need the classical trajectory is in the calculation of 9' in the Eq. ([IT)) . However, 
ip and y^{S4>'^) have the same time dependence and they only differ by a normalization which is given by Eq. ()79p . 
One can simply adopt our previous results, Eq. (|54p . for the final curvature perturbations but with \/ {6ip'^) as the 
classical trajectory and take into account the extra modifying factor ■^l/2e^ mentioned above. Using Eq. ()79|) and 
Eq. ([8T)) in Eq. ([54)) . the amplitude of curvature perturbation therefore is 



7^f ~7^^ 



1- T- 



3/2 /947r^-5/6^ 



kV'' /247re 



i, 



kcj, \ 9l3y* 



(99) 



in which we considered /3 ~ /3, fee = c^^Hq and neglected the r dependence in the exponent. Similar to the analysis 
in classical case, the induced curvature perturbations becomes 



AT'tz e^^^^ fk 



Vno cfP^'^ \k. 



(100) 



As in [2l|, [2^, Eq. ()100p shows that the induced curvature perturbations from the entropy perturbations has the 
power spectrum ex k^ and is suppressed compared to the primordial curvature perturbations by the factor e~^^''. As 
N^ <^ Nc and e~^^'' ~ 10~^* one concludes that -^^ <^ 1 for any reasonable values of the coupling g. This indicates 
that there would be no large scale curvature perturbations once the quantum back-reactions are turned on during the 
phase transition. Physically this means that the cumulative quantum back-reactions of very small scale modes, modes 
which become tachyonic during the phase transition but remained sub-horizon during entire inflation, shuts off the 
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background ip instability, forcing (f> and ^jJ to their global minima ending inflation quickly. The tachyonic instability 
of entropy perturbations are lifted and they can not produce appreciable large scale curvature perturbations during 
the phase transition. Our results can be compared with the findings of [20, |3l|, 1221 , see also [33, |3^ , although these 
works concern about the second order curvature perturbations effects. 

We can also compare our results for AV-ji with the result obtained in [44I using the 6N formalism. Noting that 

Vno - 5V47r^/3r2, we obtain AVn - £^^"^^(77)^ whereas in El] they have AVn ~ e0^^^^(^)^- The agreements 

between these two results are good (the two estimations of y'^ differ by one order of magnitude in 3Nc ~ 78 orders 
of magnitude). 

One concern may be the overproduction of primordial black holes in this model. This question was studied in [55. 
As they showed, with /3 ^ 1, the standard hybrid inflation model is safe under overproduction of primordial black 
holes. 

VI. CONCLUSIONS 

In this paper the possibility of producing large scale curvature perturbations induced from the entropy perturbations 
during the waterfall phase transition in hybrid inflation are studied. We have shown that whether or not appreciable 
amounts of large scale curvature perturbations are produced depend crucially on the competition between classical and 
quantum mechanical back-reactions to terminate inflation. If one considers only the classical back-reaction effects, 
one obtains a significant large scale curvature perturbations which completely dominate over the initial curvature 
perturbations. The induced large scale curvature perturbations would be highly blue-tilted with n-R. ~ 4 as in 
[35|. However, we have shown that the quantum- mechanical back-reactions of the waterfall field inhomogeneities 
produced during the phase transition dominate before the classical-back- reaction becomes important. In the Hartree 
approximation we found that the quantum back-reactions shuts off the classical tachyonic instability very efficiently 
terminating inflation as well as curvature perturbations evolutions quickly after phase transition. We have shown that 
the main contribution to quantum back-reactions comes from the cumulation of the very small scales inhomogeneities, 
modes which arc tachyonic during the phase transition but remain sub-horizon during entire inflationary period. We 
also made the interesting observation that in standard hybrid inflation where the waterfall field rapidly freezes to 
■0 = at the background level, one can use ip ~^ yi^W) ^^ the effective classical trajectory. In summary, the 
quantum back-reactions have two crucial effects. First they determine the end of infiation given by Eq. (|79|) and Eq. 



(fSTI) . Second, as mentioned above, they provide the effective classical trajectory via V' — > y/Wip) ■ 

Although we have presented the analysis here only for the standard hybrid inflation, but we believe that this picture 
also holds for other models of inflation where there are sharp phase transitions at the end of inflation. This includes 
models of brane infiation where infiation ends abruptly due to tachyon formation once the distance between the brane 
and anti-brane reaches a critical value. However, it would be interesting to see what happens in models of infiation, 
such as in double inflation [3a| where there is a mild phase transition in flelds evolution during early stages of inflation. 
Since the difference Nc — N^, in Eq. (llOOp is not very large, an appreciable amount of curvature perturbations can be 
created even when the quantum mechanical back-reactions are taken into account. This in turn can produce features 
in curvature perturbations such as in models [37H39J . It would be interesting to see the observational effects of the 
phase transitions during inflation as considered e.g. in (40ll41|. 
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Appendix A: Properties of Airy functions 

Here we summarize some important properties of the Airy functions which are used in the main text. 
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To the second leading order the Airy functions of second kind and their derivatives have the following asymptotic 
expansion for large arguments 

Bi(z)^^-i/2z-i/V(l + cir') (Al) 

Bi'(z)^7r-i/V/V(l + diri) (A2) 

(A3) 

in which 

C=lz"' (A4) 

and ci = 5/72 and di = —7/5 ci. By using the above equations one can estimate the derivative of Airy functions to 
the second leading order as 

Bi'(z) ^ V^ Bi(z) {l - ^^ . (A5) 

By using the above approximation for the derivative of Airy function one can find 

Bi(az)2dz ^ -l^(z)-i/2 Bi(az)2 , (A6) 

and also 

zm{azfdz~-krj:r{zY'^'Qi{azf. (A7) 

2a'^'^ 

Similarly, one can show that in the leading order the following relation also holds 

/"z"Bi(az)2dz~ ^^(z)"-i/2 Bi(az)^ (A8) 

Also one can find the following useful relation 

z" e"^Bi(az)2dz~ ^-^{zT'^^'^ Bi{azf . (A9) 

2 a"^/^ 



which was used to estimate the integral in Eq. ((52)) . 

Finally there is another simple but important relation for second derivative of Airy functions 

Bi"(z) = zBi(z) . (AlO) 



References 



[1] A. H. Guth, "The Inflationary Universe: A Possible Sohition To The Horizon And Flatness Problems," Phys. Rev. D 23, 
347 (1981); 

A. D. Linde, "A New Inflationary Universe Scenario: A Possible Solution Of The Horizon, Flatness, Homogeneity, Isotropy 
And Primordial Monopole Problems," Phys. Lett. B 108, 389 (1982); 

A. Albrecht and P. J. Steinhardt, "Cosmology For Grand Unified Theories With Radiatively Induced Symmetry Breaking," 
Phys. Rev. Lett. 48, 1220 (1982). 

[2] E. Komatsu et al., "Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpreta- 
tion," arXiv:1001.4538 [astro-ph.CO]. 

[3] A. D. Linde, "Hybrid inflation," Phys. Rev. D 49, 748 (1994) [arXiv:astro-ph/9307002]. 



18 

[4] E. J. Copeland, A. R. Liddle, D. H. Lyth, E. D. Stewart and D. Wands, "False vacuum inflation with Einstein gravity," 

Phys. Rev. D 49, 6410 (1994) [arXiv:astro-ph/9401011]. 
[5] J. Garcia-Bellido, A. D. Linde and D. Wands, "Density perturbations and black hole formation in hybrid inflation," Phys. 

Rev. D 54, 6040 (1996) [arXiv:astro-ph/9605094]. 
[6] A. Taruya and Y. Nambu, "Cosmological perturbation with two scalar flelds in reheating after inflation," Phys. Lett. B 

428, 37 (1998) [arXiv:gr-qc/9709035]. 
[7] B. A. Bassett, D. I. Kaiser and R. Maartens, "General relativistic preheating after inflation," Phys. Lett. B 455, 84 (1999) 

[arXiv:hep-ph/9808404]. 
[8] F. Finelli and R. H. Brandenberger, "Parametric ampliflcation of gravitational fluctuations during reheating," Phys. Rev. 

Lett. 82, 1362 (1999) [arXiv:hep-ph/9809490] ; F. Finelli and R. H. Brandenberger, "Parametric amplification of metric 

fluctuations during reheating in two field models," Phys. Rev. D 62, 083502 (2000) [arXiv:hep-ph/0003172]. 
[9] J. P. Zibin, R. H. Brandenberger and D. Scott, "Backreaction and the parametric resonance of cosmological fluctuations," 

Phys. Rev. D 63, 043511 (2001) [arXiv:hep-ph/0007219]. 
[10] B. A. Bassett and F. Viniegra, "Massless metric preheating," Phys. Rev. D 62, 043507 (2000) [arXiv:hep-ph/9909353]. 
[11] B. A. Bassett, C. Gordon, R. Maartens and D. I. Kaiser, "Restoring the sting to metric preheating," Phys. Rev. D 61, 

061302 (2000) [arXiv:hep-ph/9909482]. 
[12] K. Jedamzik and G. Sigl, "On metric preheating," Phys. Rev. D 61, 023519 (2000) [arXiv:hep-ph/9906287]. 
[13] A. R. Liddle, D. H. Lyth, K. A. Malik and D. Wands, "Super-horizon perturbations and preheating," Phys. Rev. D 61, 

103509 (2000) [arXiv:hep-ph/9912473]. 
[14] P. Ivanov, "On generation of metric perturbations during preheating," Phys. Rev. D 61, 023505 (2000) [arXiv:astro- 

ph/9906415]. 
[15] T. Tanaka and B. Bassett, "Application of the separate universe approach to preheating," arXiv:astro-ph/0302544. 
[16] S. Tsujikawa and B. A. Bassett, "When can preheating affect the CMB?," Phys. Lett. B 536, 9 (2002) [arXiv:astro- 

ph/0204031]. 
[17] S. Tsujikawa, D. Parkinson and B. A. Bassett, "Correlation-consistency cartography of the double inflation landscape," 

Phys. Rev. D 67, 083516 (2003) [arXiv:astro-ph/0210322]. 
[18] J. R. Bond, A. V. Frolov, Z. Huang and L. Kofman, "Non-Gaussian Spikes from Chaotic Billiards in Inflation Preheating," 

Phys. Rev. Lett. 103, 071301 (2009) [arXiv:0903.3407 [astro-ph.CO]]. 
[19] K. Kohri, D. H. Lyth and C. A. Valenzuela- Toledo, "On the generation of a non-gaussian curvature perturbation during 

preheating," JCAP 1002, 023 (2010) [arXiv:0904.0793 [hep-ph]]. 
[20] L. P. Levasseur, G. Laporte and R. Brandenberger, "Analytical Study of Mode Coupling in Hybrid Inflation," 

arXiv: 1004. 1425 [hep-th]. 
[21] D. H. Lyth, "Issues concerning the waterfall of hybrid inflation," arXiv: 1005.2461 [astro-ph.CO]. 
[22] J. Fonseca, M. Sasaki and D. Wands, "Large-scale Perturbations from the Waterfall Field in Hybrid Inflation," 

arXiv: 1005.4053 [astro-ph.CO]. 
[23] L. Randall, M. Soljacic, A. H. Guth, "Supernatural inflation: Inflation from supersymmetry with no (very) small parame- 
ters," Nucl. Phys. B472, 377-408 (1996). [hep-ph/9512439]. 
[24] E. J. Copeland, S. Pascoli and A. Rajantie, "Dynamics of tachyonic preheating after hybrid inflation," Phys. Rev. D 65, 

103517 (2002) [arXiv:hep-ph/0202031]. 
[25] C. Gordon, D. Wands, B. A. Bassett and R. Maartens, "Adiabatic and entropy perturbations from inflation," Phys. Rev. 

D 63, 023506 (2001) [arXiv:astro-ph/0009131]. 
[26] G. N. Felder, J. Garcia-Bellido, P. B. Greene, L. Kofman, A. D. Linde and I. Tkachev, "Dynamics of symmetry breaking 

and tachyonic preheating," Phys. Rev. Lett. 87, 011601 (2001) [arXiv:hep-ph/0012142]. 
[27] G. N. Felder, L. Kofman and A. D. Linde, "Tachyonic instability and dynamics of spontaneous symmetry breaking," Phys. 

Rev. D 64, 123517 (2001) [arXiv:hep-th/0106179]. 
[28] L. Kofman, A. D. Linde and A. A. Starobinsky, "Reheating after inflation," Phys. Rev. Lett. 73, 3195 (1994) [arXiv:hcp- 

th/9405187]; 

L. Kofman, A. D. Linde and A. A. Starobinsky, "Towards the theory of reheating after inflation," Phys. Rev. D 56, 3258 

(1997) [arXiv:hep-ph/9704452]. 
[29] J. H. Traschen and R. H. Brandenberger, "PARTICLE PRODUCTION DURING OUT-OF-EQUILIBRIUM PHASE 

TRANSITIONS," Phys. Rev. D 42, 2491 (1990). 
[30] Y. Shtanov, J. H. Traschen and R. H. Brandenberger, "Universe reheating after inflation," Phys. Rev. D 51, 5438 (1995) 

[arXiv:hep-ph/9407247]. 
[31] N. Barnaby and J. M. Cline, "Nongaussianity from Tachyonic Preheating in Hybrid Inflation," Phys. Rev. D 75, 086004 

(2007) [arXiv:astro-ph/0611750]; 

N. Barnaby and J. M. Cline, "Nongaussian and nonscale-invariant perturbations from tachyonic preheating in hybrid 

inflation," Phys. Rev. D 73, 106012 (2006) [arXiv:astro-ph/0601481]. 
[32] A. Mazumdar and J. Rocher, "Particle physics models of inflation and curvaton scenarios," arXiv:1001.0993 [hep-ph]. 
[33] K. Enqvist, A. Jokinen, A. Mazumdar, T. Multamaki and A. Vaihkonen, "Non-Gaussianity from Preheating," Phys. Rev. 

Lett. 94, 161301 (2005) [arXiv:astro-ph/0411394]; 

K. Enqvist, A. Jokinen, A. Mazumdar, T. Multamaki and A. Vaihkonen, "Non-gaussianity from instant and tachyonic 

preheating," JCAP 0503, 010 (2005) [arXiv:hep-ph/0501076]. 
[34] K. Enqvist and A. Vaihkonen, "Non-Gaussian perturbations in hybrid inflation," JCAP 0409, 006 (2004) [arXiv:hcp- 

ph/0405103]. 



19 

[35] J. O. Gong and M. Sasaki, "Curvature perturbation spectrum from false vacuum inflation," JCAP 0901, 001 (2009) 

[arXiv:0804.4488 [astro-ph]]. 
[36] J. Silk and M. S. Turner, "Double Inflation," Pliys. Rev. D 35, 419 (1987); 

D. Polarski and A. A. Starobinsky, "Spectra of perturbations produced by double inflation with an intermediate matter 

dominated stage," Nucl. Phys. B 385, 623 (1992); 

D. Polarski and A. A. Starobinsky, "Isocurvature perturbations in multiple inflationary models," Phys. Rev. D 50, 6123 

(1994) [arXiv:astro-ph/9404061]; 

D. Langlois, "Correlated adiabatic and isocurvature perturbations from double inflation," Phys. Rev. D 59, 123512 (1999) 

[arXiv:astro-ph/9906080]. 
[37] M. Joy, V. Sahni and A. A. Starobinsky, "A New Universal Local Feature in the Inflationary Perturbation Spectrum," 

Phys. Rev. D 77, 023514 (2008) [arXiv:071 1.1585 [astro-ph]]. 
[38] D. Battefeld, T. Battefeld and A. C. Davis, "Staggered Multi-Field Inflation," JCAP 0810, 032 (2008) [arXiv:0806.1953 

[hep-th]]. 
[39] D. Battefeld, T. Battefeld, H. Firouzjahi and N. Khosravi, "Brane Annihilations during Inflation," arXiv: 1004. 1417 [hep- 
th]. 
[40] D. Parkinson, S. Tsujikawa, B. A. Bassett and L. Amendola, "Testing for double inflation with WMAP," Phys. Rev. D 

71, 063524 (2005) [arXiv:astro-ph/0409071]. 
[41] M. Joy, A. Shafieloo, V. Sahni and A. A. Starobinsky, "Is a step in the primordial spectral index favored by CMB data 

?," JCAP 0906, 028 (2009) [arXiv;0807.3334 [astro-ph]]. 
[42] J. O. Gong and M. Sasaki, "Waterfall field in hybrid inflation and curvature perturbation," arXiv:1010.3405 [astro-ph. CO]. 
[43] D. Wands, K. A. Malik, D. H. Lyth et al., "A New approach to the evolution of cosmological perturbations on large scales," 

Phys. Rev. D62, 043527 (2000). [astro-ph/0003278] . 
[44] D. H. Lyth, "The contribution of the hybrid inflation waterfall to the primordial curvature perturbation," 

[arXiv: 1012.4617 [astro-ph.CO]]. 



